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A profile  of  the  ice  cover  in  the  southern  Beaufort  Sea 
was  obtained  h\  tin*  submarine  USS  <■!  R\AI1I)  in  April  107h,  using  a 
narrow-beam  upward-looking  sonar.  The  1 ri()0  km  profile  consisted 
of  three  legs,  of  which  the  long  south-north  and  east-west  legs 
intersected  near  the  Caribou  camp  of  the  AlItlKX  experiment,  A 
statistical  analysis  was  rallied  out  over  contiguous  'jO  km 
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sections  to  yield  probability  density  functions  of  the  drafts  of 
ice  and  of  level  ice,  distributions  of  keel  spacing*  and  drafts, 
and  the  frequencies  and  widths  of  leads.  Two  distinct  types  of 
ice  cover  were  found.  The  first,  nearest  the  coast  in  the  south 
and  west  of  the  experimental  area,  consisted  of  heavily  ridged  ice 
with  mean  drafts  of  up  to  5.1  m.  The  rest  of  the  track  (1200  km) 
consisted  of  a homogeneous  ice  cover  with  a mean  draft  of  3.7  m« 
The  percentage  of  thin  ice  varied  greatly  between  sections,  with 
a range  of  0.4  to  12.3^  for  ice  of  0-1  m draft.  J.evel  ice,  with 
a local  gradient  of  less  than  1 in  40,  made  up  5b^n  of  the 
homogeneous  cover,  with  a preferred  draft  of  2.7  to  2.9  m.  Keel 
spacing*  obeyed  a negative  exponential  distribution,  with  a 
deficit  at  small  spacings  due  $.0  a keel  shadowing  effect  and  a 
surfeit  at  very  large  spacings  due  to  the  contribution  of 
polynyas.  The  draft  distribution  of  keels  was  also  a negative 
exponential,  differing  from  the  distribution  of  Hibler  cJ  a_j_ 
(1972),  probably  because  the  narrow  beam  records  a complex 
structure  for  every  keel.  The  coastal  ice  had  a greater 
frequency  and  mean  draft  of  keels.  Maximum  keel  draft  was 
31.12  m.  The  mean  separation  of  lends  was  212  m,  with  almost  all 
lends  being  less  than  50  m in  width. 
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SUMMARY 

A profile  of  the  ice  cover  in  the  southern  Beaufort  Sea  was 

obtained  by  the  nuclear  submarine  USS  GURNARD  from  7 to  10  April 

197b, 

using  a narrow-beam  upward-looking  sonur.  The  1400  km  profile 

consisted  of  three  legs,  of  which  the  long  south-north  and  east-west 

legs 

intersected  near  the  Caribou  ramp  of  the  AIDJEX  experiment. 

After  initial  processing  at  the  Arctic  Submarine  Laboratory,  San  Diego, 

the 

duta  were  passed  vin  the  AIIXJKX  Office  to  Scott  Polar  Research 

Institute  for  statistical  analysis. 

The  analysis  was  carried  out  over  contiguous  10  km  sections 

of 

the 

profile,  and  yielded  the  following  significant  results:- 

1 

i I 

J 


I.  I'rohubi  1 i t \ density  function  of  ire  draft.  The  profile  was  found  to 
consist  of  Ihret*  distinct  types  of  ice  cover.  Type  1,  nearest  the  coast 
and  extending  from  the  100  m isobath  northwards  to  71°A0'N  at  IVi°1VW, 
consisted  of  heavily  ridged  ice  with  a mean  draft  of  ?».2-A.I  m.  Type  'J 
comprised  the  major  part  of  the  profile  (1200  km),  mid  was  shown 
statistically  to  be  a homogeneous  ice  cover  with  the  following  distribution 
of  ice  drafts:-  0-0. 1 m,  1%;  0.A-2  m,  1 1>* »■;  2- A m,  over  r>  tu,  1 7** <•; 

mean  draft  A.h7  m.  Type  A was  another  short  stretch  of  track  at  the  far 
western  end  of  the  profile  (1A1°W  to  IA4°1VW  at  72<>,i2,N)  which  again 
consisted  of  more  heavily  ridged  ice  of  mean  draft  'i.A-A.h  in.  The 
proportion  of  thin  ice  was  the  most  variable  quantity  in  the  distributions, 
having  a range  0.1 -A. AY  (0-0.r>  m draft)  and  0.1|-12.AA  ( 0—  l m draft).  The 
cumulative  probability  distributions  for  two  200  km  sections  crossing  the 
A I HI KX  comps  were  calculated  for  application  to  the  A11VIKX  model. 


2.  level  ire.  Ice  of  less  than  1 in  hQ  local  gradient  - i.e.  mainly 
undeformed  floes  eont ami na t ed  with  a lew  keel  bottoms  - was  found  to  comprise 
Ati‘<  of  the  ice  in  the  type  2 cover.  The  most  probable  draft  of  level  ice 
was  2. 7-2.8  m,  with  secondary  peaks  at  2.0-2. 1 m and  A.1-A.2  m.  Thin  ice 
had  a peak  at  0. 8-0.0  m. 

A.  Keel  spacing*.  l’he  distribution  of  keel  spocings  was  found  to  obey  a 
negative  exponential  distribution,  with  departures  at  small  spaeings  (a 
deficit  line  to  a "keel  shadowing"  effect,  where  shallow  keels  are  concealed 
by  neighbouring  deep  keels)  and  at  large  spaeings  exceeding  A00  m (a  surfeit 
due  to  the  presence  of  leads  and  polynyas). 

’i  ■ Keel  drafts.  fhe  distribution  of  drafts  did  not  obey  the  theory  of 
llibler  e_t_  id  {To 72 ) which  predicts  l‘(h)  exp(-bh“’),  hut  rather  a simpler 

distribution  of  form 

l’(h)  dh  = H exp(-bh)  dh 

with  11,  b as  parameters.  Hitherto  only  sails  have  been  found  to  obey 
this  relationship,  and  it  is  postulated  that  the  cause  is  the  very  narrow 
be<unwidth  of  the  transducer,  which  enables  it  to  probe  into  the  fine 
structure  of  a keel  mid  thus  perceive  multipit*  secondary  "keels"  where  a 
wide-beoin  sounder  sees  only  a single  entity.  Again  types  1 and  A have  heavier 
ridging  than  type  2.  For  keels  deeper  than  0 m,  type  2 ice  has  a mean  of 
1 .AO  keels  per  km  with  a mean  draft  of  11. 7*^  ni ; types  1 mid  A have  up  to 
A.l  keels  per  km  with  a me  mi  draft  of  12.1-12.7  m.  Maximum  keel  draft  was 
A 1.12  in,  with  only  one  keel  over  At)  m in  draft. 

A.  Leads  and  polynyas.  A lead  was  defined  as  a continuous  sequence  of 
depth  points  where  no  point  exceeds  l m in  draft.  The  average  spacing  of 
leads  was  212  m over  the  whole  profile,  but  almost  all  leads  were  less  than 
A0  m in  width.  The  mean  spacing  of  leads  over  A0  m wide  was  10. A kui,  and 
of  leads  over  A00  m wide  was  2A7  km. 
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During  the  period  7 to  10  April  l')7b  the  line  lour  submarine 
1'SS  GPRNARD  (SSN-tit)2)  obtained  a sonar  profile  of  length  1A00  km 
under  the  Beaufort  Sea  ice  cover  in  the  vicinity  of  the  A1RJEX 
(Arctic  Ice  Dynamics  Joint  Experiment')  mam  camp.  Figure  1 shows 
the  route  followed  b\  GPRNARD,  comprising  south-north  (Ol’Q)  ami 
east-west  (Ill’S)  legs  intersecting  near  the  Caribou  c.unp  (C), 
together  with  a connecting  leg  yil . (iCllNARD  was  equipped  with  a 
high-frequency , nar  row-beam,  upward-looking  sonar  installed  by  the 
Arctic  Submarine  laboratory.  Naval  Undersea  Center,  Sun  Diego.  The 
sonar  fed  its  output  into  a signal  processing  system  that  digitised 
the  range  to  the  ice  underside,  subtracted  tins  from  tile  transducer 
depth  and  thus  generated  a digital  magnetic  tape  of  ice  drafts  with 
a nominal  resolution  of  0.03  m.  Mechanical  limitations  within  the 
system  reduced  the  absolute  accuraci  of  ice  draft  to  + 0.3  m,  with 
a demonstrated  standard  deviation  for  smooth  ice  of  0.00  m.  Other 
details  of  the  sonar  system  and  of  the  submarine's  depth  and  speed 
remain  classified.  The  Arctic  Submarine  Laboratory  did,  however, 
supply  information  on  the  "surface  beam  diameter"  (diameter  of  spread 
of  the  sonar  beam  at  the  surface,  a function  of  hcumwidth  ;uul  cruising 
depth)  and  the  "ping  spacing"  (bor i /ontal  distance  between  successive 
sound  pulses,  a function  of  ping  frequency  and  submarine  speed). 

Normal  1\  the  pmg  spacing  was  between  1.3  and  1.3  m,  and  the  surface 
beam  diameter  over  almost  the  whole  track  was  3.17  m,  implying  a very 
narrow  beam  of  less  than  3 width. 

Initial  processing  of  the  tapes  was  done  at  the  Arctic  Submarine 
Laboratory  using  a Pnivac  1108-1110  computer.  Corrected  depth  duta 
were  merged  with  positional  information  to  give  a final  tape  in  BCD 
format  which  was  forwarded  to  the  AIDJKX  Project  Office  in  Seattle  for 
analysis.  In  turn  the  A1IXJEX  Office  forwarded  the  tape  to  the  Scott. 
Polar  Research  Institute  so  that  the  data  could  be  analysed  using  the 
same  criteria  and  definitions  as  those  employed  in  the  analysis  of  data 
from  IIMS  SOVEREIGN  (Wadbams  1 ‘>77 »» , b ; Wadhams  and  Lowrv,  1077)-  This 
report  gives  the  results  of  such  an  analysis. 


J.  DATA  PROCESSING 


The  initial  processing  at  the  Arctic  Submarine  Laboratory 
deleted  spurious  profile  points  caused  by  multiple  echoes,  fish,  air 
bubbles  etc.,  and  set  these  to  -10.0  ft.  Isolated  spurious  points 
were  then  regenerated  by  linear  interpolation.  Appendix  A describes 
this  stage  of  processing  in  more  detail,  including  the  criteria  for 
identifying  spurious  points.  The  data  file  created  at  San  Diego 
consisted  of  a series  of  "blocks"  each  containing  about  (>0  data  points. 
The  blocks  were  separated  b\  single  lines  which  usually  contained  all 
zeroes,  but  which  at  periodic  intervals  contained  a position  fix 
(i.e.  latitude  and  longitude  values).  These  position  fixes  thus  split 
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FIGURE  1 (above).  Route  of 
USS  GURNARD,  7-10  April  107b 


FIGURE  2 (left).  Division  of 
track  into  'ill  km  sections. 
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tin*  data  file  into  "intervals"  for  which  the  distance  travelled  is 
known . 

The  analysis  at  Scott  Polar  Research  Institute  was  carried 
out  on  the  I'niversitv  of  Cambridge  IBM  T7t>/ 1 computer  using  an 
Algol  b8C  program  which  is  listed  and  described  in  Appendix  H.  The 
program  split  the  input  file  into  "sections".  Each  section  contained 
sufficient  intervals  to  make  up  ”>0  km  of  data,  and  statistics  were 
computed  for  each  of  the  sections.  figure  2 shows  the  27  sections 
involved  (with  an  18  km  end-of-file  gap  between  sections  12  mid  17)  and 
their  precise  positions  and  lengths  are  given  in  table  1. 


fable  1.  Positions  of  the  70  km  sections. 


Section  latitude  Longitude  True  length 

km 


1 

Start 

s 70 

' 7 V 

1 1 . 2" 

144° 

17' 

14.4" 

72.07 

Ends 

71 

01 

'»0.2 

144 

17 

18.0 

• > 

Ends 

71 

10 

18 . 2 

144 

17 

47.0 

70.71 

A 

Ends 

71 

18 

01.2 

144 

14 

20 . 4 

70.72 

4 

Ends 

21 

2b . 4 

144 

14 

78.8 

71.01 

7 

End  s 

-•) 

12 

• >•>  •> 

144 

17 

72 . 8 

12.17 

t> 

f Jills 

71 

20 

07.2 

144 

18 

Ob  .0 

71.02 

f 

Ends 

71 

48 

18.0 

144 

20 

02 . 4 

72.47 

8 

f inds 

7-i 

!(■ 

Vi.2 

144 

*>•» 

70. t) 

72.01 

0 

Ends 

7'* 

•’ll 

7b . k 

144 

')•) 

27.8 

70.02 

10 

Ends 

71 

08 

01.  '1 

144 

18.0 

70.77 

1 1 

Ends 

71 

11 

12 . 2 

141 

47 

20.4 

71 . 4 b 

12 

Ends 

7'* 

12 

11.4 

1 42 

48 

21. b 

72 . 78 

n 

Start 

s 7** 

10 

28.8 

142 

lb 

70 . 4 

72.40 

Ends 

7'# 

11 

'il.  2 

141 

20 

22.8 

14 

Ends 

71 

12 

00.  b 

1 40 

21 

20.4 

12.04 

ii 

Ends 

71 

28 

'1 5 • 0 

110 

72 

24.0 

71.07 

lb 

finds 

71 

04 

0 '» . 2 

118 

47 

47.4 

72.80 

17 

Ends 

•70 

’ll) 

10.2 

178 

11 

01.8 

71.70 

18 

finds 

’•0 

'4  2 . b 

110 

40 

40.8 

72.72 

10 

Ends 

•’ll 

28.8 

141 

21 

70.  b 

77.08 

20 

finds 

*?•> 

k2 

40.2 

142 

1b 

00 . 0 

70.04 

21 

finds 

*70 

k5 

10.2 

144 

2t> 

77.  b 

70 . 4 1 

22 

Star  t 

^ *70 

47 

22 . 0 

144 

48 

71.2 

72.08 

finds 

47 

1 4 . 4 

1 4 b 

27 

11  .b 

21 

Ends 

-»o 

42 

48.0 

147 

70 

18.0 

72.00 

2k 

finds 

7* » 

42 

01.4 

140 

72 

17.2 

71.7b 

25 

finds 

’•1 

21.8 

111 

02 

72.4 

70.00 

2b 

finds 

't) 

k2 

01.2 

112 

78 

48.  b 

17 . 20 

27 

finds 

70 

'il 

2° . ■* 

1 14 

14 

72 . 4 

72.08 
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Eoi 
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th*'  A 1 n «» 1 program  therefore  h«*l  to  calculate  tin  "interpolation  length", 
the  trur  spacing  between  depth  points,  l»\  dividing  the  length  of  the 
interval  (calculated  assuming  a (treat  Circle  track  between  the  position 
fixes  for  the  beginning  and  end  of  the  interval)  by  the  number  of  points 
in  the  interval.  The  contribution  made  by  each  interval  to  the  overall 
statistics  for  a 'll)  km  section  was  then  always  weighted  In  the  inter- 
polation length,  so  that  tin*  resulting  statistics  are  unbiased  with 
respect  to  horizontal  length. 

The  statistics  generated  are  described  and  interpreted  m the 
remainder  of  this  report,  and  the  full  numerical  values  are  given  in 
Appendix  (’.  The  program  produced  two  types  of  output  f i 1 e : — 

(a)  a line  printer  listing  of  the  statistics,  with  details  of 
the  intervals  making  up  the  section  (Appendix  C); 

(b)  a file  containing  three  frequency  tables  for  each  section 
(ice  draft  distribution,  level  ice  and  ridges)  from  which 
histograms  were  produced  using  a separate  program.  These 
appear  in  the  body  of  the  report. 

figure  3 shows  a small  part  of  the  depth  data  plotted  out. 

The  two  most  obvious  features  are  a high-frequency  noise  superimposed 
on  the  supposedly  smooth  ice  bottom  contour,  and  the  occasional 
shallow  depth  point  occurring  within  the  structure  of  a pressure  ridge. 
The  noise  is  undoubtedly  a feature  of  the  recording  system  and,  since 
it  is  random,  it  does  not  have  a serious  effect  on  probability 
densities  of  draft,  although  it  will  produce  tin  anomalously  high  number 
of  very  small  "pressure  ridges"  in  the  statistics.  The  shallow  depth 
points  are  probably  real  features,  caused  by  the  very  narrow  sonar 
beam  probing  into  fissures  and  crannies  within  the  loose  block  structure 
of  the  ridge.  The  effect  of  these  points  may  be  quite  large  in 
causing  a single  pressure  ridge  to  appear  as  multiple  ridges  in  the 
statistics. 


>.  PHI  INABILITY  DENSITY  OF  ICE  DRAFT 


3.1.  Def  ini  t i».:i 


The  probability  density  function  l’(h)  of  draft  h is  defined  such 
that  P(h)  dh  is  the  probability  ihat  a random  point  oil  the  ice  underside 
has  u draft  between  li  and  (h  + dh).  P(h)  should  really  be  expressed  in 
the  form  P(h,  t)  since  it  is  a function  of  time  as  well  as  of  position. 
Further,  although  P can  be  stochastically  defined  at  a point  an 
operational  definition  requires  a profile  to  be  taken  over  a finite  length 
scale  in  order  to  arrive  at  an  unbiased  estimate  of  P.  This  length  scale 
must  be  large  enough  to  give  a good  estimator  of  P while  small  enough  for 
the  distribution  not  to  change  significantly  within  its  compass.  We  have 
chosen  *>0  km,  blit  some  analyses  have  been  done  over  shorter  (17  km)  and 
longer  (200  km)  length  scales,  where  necessary. 

P(h)  is  related  via  the  mean  density  of  the  ice  to  the  thickness 
probability  density  function  g(h)  of  Thorndike  al.  (1973).  g(h)  is 
important  as  ,ui  input  parameter  to  various  models  of  Arctic  Ocean  ice 


B 


dynamics  arid  thermodynamics.  These  include  the  A1DJEX  model 
(Coon  et  al,  197**)*  in  which  an  initial  ice  thickness  distribution 
develops  by  thermodynamic  growth  and  decuy  and  is  continuously 
redistributed  b\  pressure  ridge  building  and  wind-driven  divergence; 
mid  the  viscous-plastic  continuum  model  of  Hibler  (1977)*  in  which 
a pressure  term  is  parameter i sed  using  the  mean  ice  thickness  and 
percentage  ice  cover.  The  most  sensitive  part  of  l’(h)  is  the  thin 
ice  component,  since  it  has  been  shown  (badgley,  1996;  Maykut,  197h ) 
that  most  of  the  heat  flow  from  ocean  to  atmosphere  in  the  Arctic 
occurs  through  ice  of  draft  less  than  1 m;  this  is  also  the  ice 
component  which  is  most  readily  available  for  ridge  building. 


"5.2.  Results 


Figure  shows  P(h)  plotted  for  all  27  sections  of  the  50  km 
length  scale,  using  a depth  increment  of  10  cm;  the  numerical  data 
on  which  these  plots  are  based  are  given  in  Appendix  C.  The  general 
nature  of  all  the  plots  is  similar:-  an  initial  peak,  due  to  thin 
ice  in  leads  and  polynyas;  a second,  broader  peak  due  mainly  to 
undeformed  first-  and  multi-year  ice;  and  a tail  due  to  ice  in 
ridges  and  hummocks.  There  is  some  variation  from  section  to 
section,  especially  in  the  extent  of  thin  ice  present. 

To  display  these  variations  more  clearl\  P(h)  was  integrated 
over  four  depth  intervals,  which  can  be  loosely  defined  as  "thin 
ice"  (0-0.5  m) ; "young  ice"  (0.5  -2m);  "level  ice"  (2-5  m)  and 
"ridged  ice"  (greater  than  5 m) . The  separation  of  types  is  not 
perfect  - parts  of  ridges,  for  instance,  may  appear  in  the  "level 
ice"  category  - but  the  categories  are  indicative  of  changes  in  the 
nature  of  the  ice  cover.  The  results  are  given  in  Table  II.  The 
intervals  were  chosen  so  as  to  give  a direct  comparison  with  the 
data  of  Wadhams  ( 1977b)  from  the  heavily  ridged  offshore  zone  to 
the  north  of  Greenland  and  Ellesmere  Island.  Wadhams  found  that 

the  "thin  ice  peak"  in  his  probability  density  functions  usually 
occurred  at  less  than  0.5  m draft,  hence  his  choice  of  intervals, 
but  the  present  results  (fig.  A)  usually  show  the  peak  at  between 
0.5  and  1 m,  presumably  because  the  profiles  were  done  later  in  the 
winter  (April  compared  to  October  for  SOVEREIGN)  so  that  the  ice  in 
polynyas  is,  on  average,  thicker.  Thu3  we  have  also  added  a 0-1  tn 
category  in  Table  II  to  include  all  of  the  polynya  ice. 

The  results  show  a remarkable  consistency  of  ice  conditions 
over  most  of  the  experimental  area.  The  exceptions  are:- 

(i)  the  percentage  of  thin  ice,  which  varies  over  a wide  range 
(0.4  to  12.35*  for  the  0-1  m band)  and  with  no  apparent 
consistency  of  trend.  The  cause  is  partly  statistical  - 
thin  ice  is  contained  in  a limited  number  of  polynyas  which 
are  distributed  non-unif ormly  along  the  submarine  track  - 
and  partly  real  in  the  sense  that  thin  ice  has  a transient 
existence  and  is  constantly  being  destroyed  by  ridge- 
building so  that  changes  in  the  wind  field  during  the  3 days 
of  the  experiment  may  cause  the  thin  ice  to  be  radically 
redistributed. 
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FIGURE  ;i . Probability 
density  function  of  ice 
draft  for  ^0— kin  sections. 
Depth  increment  10  cm. 
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Table  II 


Perron tages  of  ire  rover  i it  different  range*  of  draft 


(a)  30  kin  see  t ions 


Ser  t i on 

0-0.3 

0.3-2 

2-3 

3 

m-» 

0-1 

Mean 
d i a f t 

in 

m 

m 

m 

1 

1.2 

11 

47 

41 

4.3 

3.09 

•> 

0.4 

13 

t>0 

23 

0.8 

4 . 22 

3 

0.3 

18 

b2 

20 

4.3 

3.77 

4 

0. 1 

•> 

78 

20 

0. 4 

'! . 1 4 

3 

1.3 

/ 

71 

21 

2.8 

3 . 92 

b 

0.7 

13 

73 

13 

7-8 

3.43 

i 

0.9 

23 

<»3 

13 

9.3 

3.19 

8 

0.2 

10 

78 

12 

2.0 

3.57 

9 

0.  1 

12 

71 

17 

0.3 

3.dl 

10 

1.2 

11 

74 

14 

.»  •» 

3.47 

11 

2.3 

11 

73 

14 

4.8 

3 . 40 

12 

1.2 

9 

70 

i — 

18 

1.5 

3.74 

13 

0.3 

9 

74 

17 

0.8 

3.70 

U» 

0.2 

3 

78 

lb 

0.7 

3.78 

13 

0.3 

9 

78 

13 

0.8 

3.31 

lb 

O.b 

9 

79 

12 

1 .0 

3.38 

17 

0.8 

8 

72 

19 

3.  1 

3.89 

18 

1.8 

10 

bb 

22 

3.9 

5.87 

19 

3.3 

10 

72 

13 

3.2 

3 . b2 

20 

2.7 

8 

70 

i - 

18 

4.4 

3. 69 

21 

2.7 

17 

b4 

lb 

8.2 

3.47 

22 

0.3 

11 

73 

14 

2.b 

3.33 

23 

0.3 

4 

7b 

20 

1.1 

4.14 

24 

0.3 

18 

b3 

lb 

12.3 

5.b3 

23 

0.3 

4 

73 

21 

1.7 

4.17 

2b 

O.b 

8 

b3 

28 

3.0 

4 . 30 

27 

0.8 

3 

b8 

27 

2.3 

4 . b l 

Mean 

0.9 

10 

71 

19 

3.4 

5.81 

SOVEREIGN 

4.8 

4 

28 

b3 

b.  3 

7 . 22 

DREADNOUGHT 

O 

- . t 

9 

33 

33 

3.3 

4.29 

17  km  sections 

in  offshore 

zone 

1.1 

0.2 

9 

42 

48 

O.b 

3 . 38 

1.2 

3.7 

14 

47 

33 

13.3 

4.43 

1.3 

0.0 

9 

33 

38 

0. 1 

3.18 

2.1 

0.3 

13 

3b 

31 

0.9 

4 . 3 1 

1 9 

i>  I k 

0.3 

18 

37 

24 

O.b 

4.14 

2.3 

0.7 

13 

bb 

21 

0.8 

4.03 

200  km  sections 

A 

O.b 

13 

3b 

29 

3.2 

4 . 3b 

B (Caribou  n-s) 

0.8 

11 

71 

17 

3.1 

3.bH 

C 

1.0 

11 

74 

14 

2.4 

3 . 4b 

D 

O.b 

8 

7b 

lb 

1.0 

3. 98 

E 

1.7 

9 

72 

17 

3.3 

3.99 

F (Caribou  e-w) 

l.b 

10 

72 

17 

4.1 

3.71 

G 

O.b 

9 

b8 

23 

4.9 

4 . 23 
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(ii)  Ridged  ice  at  the  southernmost  (and,  to  a leaser  extent,  the 
westernmost)  end  of  the  profile  is  significantly  greater  in 
quantity.  The  percentage  of  ridged  ice  is  very  high  in 
section  1 and  diminishes  to  u fairly  steady  "equilibrium 
value"  by  section  3.  Clearly  the  first  2 sections  represent 
the  "offshore  province"  of  Weeks  et,  a_l  ( 1 97 1 ) * a heavily 
ridged  coastal  province  where  the  mean  onshore  tendency  of 
ice  drift  leads  to  net  convergence  and  ridge  building.  To 
a lesser  extent  sections  26-27  mark  the  outer  edges  of  this 
province  further  to  the  west.  Figure  1 shows  that  sections  1 
and  2 occur  over  the  continental  slope  (the  profile  itself 
commences  ut  the  100  m depth  contour)  and  that  section  27 
ends  just  as  this  slope  is  again  approached  off  Point  Harrow. 

The  remainder  of  the  sections  are  samples  of  what  seems  a 
very  homogeneous  ice  cover. 

The  mean  values  over  all  27  sections  are  compared  in  table  II 
with  the  mean  results  from  SOVEREIGN  (a  1000  km  profile  from  81°N  0°W 
to  8'i°50'N  70°W)  and  DREADNOUGHT  (a  560  km  profile  from  85°  to  90°N 
at  b°E  in  the  ice  of  the  Trans  Polar  Drift  Stream;  Williams  ej_  al , 
1975).  The  best  agreement  is  with  the  DREADNOUGHT  data,  although  the 
GURNARD  data  show  a somewhat  lower  mean  draft  which  can  be  ascribed  to 
the  beamwidth  of  the  DREADNOUGHT'S  echo  sounder.  Clearly  the  ice 
encountered  by  SOVEREIGN  was  far  more  heavily  ridged  and  thicker  than 
even  the  heaviest  section  of  GURNARD. 

The  50  km  gauge  is  short  enough  to  resolve  most  real  variations, 
but  to  investigate  the  rapidly  changing  ice  conditions  at  the  beginning 
of  the  profile  (sections  1 and  2)  a 17  km  gauge  was  used,  i.e.  each 
section  was  split  into  three.  The  results,  in  table  11(b),  show  a 
steady  decrease  in  the  percentage  of  ridged  ice  as  the  submarine  travels 
north  away  from  the  Alaskan  coast.  Note  the  isolated  value  of  15*3% 
for  0-1  m ice  in  section  1.2. 

Finally,  to  obtain  very  reliable  statistics  for  large  tracts  of 
the  ice  cover,  the  50  km  sections  were  combined  into  a 200-km  length 
gauge  as  shown  in  figure  2,  lettered  A to  G (A  is  150  km  only).  B 
and  F are  now  the  appropriate  sections  for  the  crossings  under  Caribou 
camp,  and  have  the  length  scale  recommended  in  the  AIDJEX.  model  and  by 
Thorndike  £t  al_  (1975)-  The  results  for  B and  F agree  very  closely. 
Again,  the  200  km  statistics  show  that  the  character  of  the  ice  is 
essentially  constnnt  over  most  of  the  track  (B  to  F),  with  an  increase 
in  moon  draft  and  percentage  of  ridged  ice  at  the  western  end  (G)  and, 
particularly,  the  southern  end  (A) . 


3.3  Statistical  reliability 

It  is  clear  from  table  11(b)  that  there  are  progressive  changes  in 
the  nature  of  the  ice  cover  over  sections  1 and  2,  and  that  these  sections 
(and  probably  26  and  27)  differ  in  nature  from  sections  3-25.  The 
question  that  remains  is  whether  the  variations  between  sections  in  3-25 
are  statistical  artefacts,  i.e.  due  to  finite  sampling  length,  or 
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whether  they  are  due  to  reul  variations,  albeit  minor,  in  the  nuture 
of  the  iee  cover.  Our  null  hypothesis  is  thut  the  ice  cover  over  a 
substantial  part  of  the  southern  Beaufort  Sea  (the  urea  sampled  by 
sections  3-25)  is  a homogeneous  cover  in  which  any  given  statistical 
parameter  tends  to  the  same  value  everywhere  if  sampled  over  a 
sufficient  length  of  track. 

The  question  of  "sufficient  length"is  a crucial  one.  For  a 
given  sampling  length,  some  parameters  are  estimated  more  accurately 
than  others.  For  instance,  50  km  of  track  usually  contains  enough 
ice  types  to  make  it  a good  estimator  of  mean  draft,  and  most  probably 
of  percentage  ridged  ice.  It  may  not  be  a good  estimator,  however, 
of  thin  ice  percentage  since  this  ice  is  contained  in  a small  number 
of  polynvas  which  may  not  happen  to  fall  uniformly  within  the  length 
gauge.  A longer  gauge  may  be  required  for  a good  estimate  of  this 
parameter,  and  also  of  such  parameters  as  mean  polvnya  spacing  or  the 
frequency  of  very  deep  pressure  ridges. 


We  can  test  for  the  homogeneity  of  the  data  using  a noil-parametric 
run  test  (Bendat  and  Piersol,  1971).  The  data  is  divided  into  n 
sections  for  which  a given  statistic  S takes  values  Sj  (j=l,n).  The 
mean  value  of  S is  calculated,  and  each  section  is  classified  as  (+)  or 
(-)  according  as  Sj  >■  S or  Sj  < S.  The  number  of  runs  of  consecutive 
(+)  or  (-)  classes  in  the  n sections  is  found  and  tested  for  significance. 
An  exceptionally  small  number  of  runs  signifies  a trend  or  a clustering 
in  S;  an  exceptionally  large  number  signifies  a factor  tending  to  cause 
alternation  of  high  and  low  S values.  The  23  50-lan  sections  (3-25)  were 
tested  in  this  way,  and  an  additional  test  was  afforded  by  30  17-km 
sections  corresponding  to  sections  3-12  (sections  1-12  were  analysed  at 
17-km  gauge  for  use  in  table  II).  The  mean  values  of  the  statistics 
of  table  II  over  3-25  and  the  results  of  the  test  are  as  follows 
(*  signifies  non-significance):- 


Mean  and  Test  - 50  km  Test  - 17  km 

standard  error  n=23  n=30 


0-0.5  m percentage  cover  1.0 
0.5-2  m " " 10.3 

2-5  m " " 72 

over  5 m " " lb. 6 

0-  m " " 3.5 

mean  draft  " " 3 .07 


+ 0.2 

2.5% 

* 

+ 0.9 

* 

* 

+ 1 

* 

* 

+ 0.7 

* 

* 

+ 0.7 

* 

* 

+ 0.0b 

* 

* 

Only  one  statistic  - the  0-0.5  m percentage  cover  - was  rejected  by 
this  test,  at  the  25%  significance  level  and  for  the  50  km  sections. 

This  implies  that  the  thin  ice  percentage  does  not  come  from  a homogeneous 
population,  i.e.  that  there  are  significant  trends  or  clusterings  in  this 
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statistic  (which  can  be  seen  bv  inspection  of  tuble  ll)  suggestive  of 
a process  acting  with  a wavelength  much  greater  than  ‘>0  km  and 
tending  to  cause  large  amounts  of  thin  ice  over  a number  of  sections 
followed  bv  small  amounts.  This  process  must  be  the  wind  stress 
field  which  causes  divergence  in  one  zone  of  the  ice  cover  and 
convergence  in  another,  on  a length  scale  of  hundreds  of  km. 

Otherwise  we  can  accept  the  hypothesis  that  the  ice  draft 
distributions  over  1150  kin  of  track  in  the  southern  Beaufort  Sea 
(5-25)  come  from  a homogeneous  ice  cover  with  constant  statistical 
properties.  The  best  values  to  take  for  the  mean  draft  and  percen- 
tages of  ice  in  various  depth  ranges  are  given  above,  together  with 
the  standard  error.  It  can  be  seen  that  this  standard  error  is 
virtually  constant  for  each  class  (except  0-0.5  m) , implying  a greater 
fractional  error  in  tin-  estimates  of  uncommon  ice  types  than  in  those 
of  common  types  - a result  to  be  expected. 

5.4  Cumulative  probability 

The  cumulative  probability  U(h)  is  defined  by 

0(h)  = /h  P(h)  dh 
o 

It  is  used  as  a major  parameter  in  the  AIDJEX  model  (Coon  et_  al_,  197*0 
where  it  is  known  as  the  "ice  thickness  distribution".  We  have 
computed  0(h)  for  the  two  200  km  sections  B and  F which  cross  the 
Caribou  camp  in  the  south-north  and  east-west  directions.  The  results 
ure  shown  in  figure  5. 

The  two  distributions  differ  slightly  but  are  similar  in  general 
shape.  The  median  depth  (0(h)  = 0.5)  is  reached  at  5.2  m and  the 

graph  is  plotted  as  far  as  0(h)  = 0.99,  which  is  reached  at  12.2  m. 

G(h)  = 0.099  is  reached  at  lb. 4 m. 


4.  LEVEL  ICE 

4 . 1 I>e  f i n i t i on 

Williams  e_t  al_  (1975)>  in  their  analysis  of  the  DREADNOUGHT  data, 
sought  a way  of  determining  the  preferred  thickness  or  thicknesses  of 
undeformed  floes,  and  the  percentage  of  the  ice  cover  occupied  by 
floes  of  this  type.  By  trial  and  error  they  decided  that  the  best 
working  definition  of  "level  ice"  is  that  the  draft  point  concerned 
should  have  draft  points  4 m to  each  side  of  it  differing  in  depth  by 
less  than  20  cm,  i.e.  a local  gradient  of  less  than  1 in  40  measured 
on  an  8 m gauge  length  across  the  point.  On  account  of  the  high- 
frequency  structure  in  the  GURNARD  profile  (fig.  5),  we  have  relaxed 
this  definition  slightly  and  we  define  a level  ice  point  as  one  whose 
draft  differs  from  a point  10  m away  to  either  side  by  less  than  25  cm, 
i.e.  a 1 in  40  gradient  in  one  direction  moving  away  from  the  point. 
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FIGURE  5 (above).  Cumulative 
probability  distribution 
G(h)  of  ice  draft  for  the 
two  200  km  sections  which 
straddle  the  Caribou  camp. 


FIGURE  7 (left).  Probability 
density  function  P(h)  of 
level  ice  draft  for  whole 
submarine  track.  Bin  size 
0.1  m.  Dotted  lines  show 
bins  between  2 and  k m draft 
after  normalisation  (see 


•'•.2  itesults 


Figure  ()  shows  probability  density  flint- 1 i «ns  of  level  ice  draft 
for  all  27  soi  l ions  at  50  km  gauge,  again  using  a 10  cm  depth  increment. 
Am  "level  ice"  occurring  beyond  about  5 m draft  must  be  on  the  bottoms 
of  ridges  and  hummocks.  The  inferences  to  be  drawn  from  these  plots 
concern: - 

(i)  positions  of  peak  or  peaks,  indicating  preferred  drafts 
of  level  ice  of  varying  ages  (young,  first-year,  second- 
year  , etc.); 

(ii)  mean  draft  of  level  ice; 

(111)  percentage  of  total  ice  cover  occupied  by  "level  ice". 

These  results  are  gathered  in  table  111. 

fhe re  are  two  methods  of  estimating  (i).  tine  is  to  plot  the 
probability  density  function  of  level  ice  over  the  whole  track  length 
(figure  7):  the  other  is  to  use  a ranking  technique  on  the  individual 
50  km  plots  to  detect  preferred  drafts.  In  table  111  this  has  been 
done  by  listing  in  order  of  sue  every  bin  with  a probability  density 
of  greater  than  0.5  plus  the  thin  ice  peak  if  it  exceeds  0.1. 

It  is  tempting  to  ascribe  the  preferred  depths  of  table  111  to 
ice  ol  varying  ages.  Clearly  the  2 . 7-2 . 8 m category  is  the  most 
popular  depth  (shown  also  by  figure  7)  with  2. 5-2. 5 m and  2.0-2. 1 m 
as  lesser  peaks.  Thick  level  ice  of  5.1  depth  and  above  is  confined 
to  regions  in  the  far  west  and  far  south  of  the  profile,  i.e.  the 
regions  of  heaviest  ridging.  However,  figure  7,  which  at  first 
sight  confirms  the  hypothesis  of  a discrete  number  of  preferred  drafts, 
in  fact  required  modification.  The  10  cm  increments  from  which 
figure  7 is  constructed  are  not  all  present  with  equal  probability, 
since  the  original  draft  data  from  San  Diego  were  quoted  in  integral 
numbers  of  (loci  feet  (0.05048  m).  Thus  some  of  the  10  cm  increments 
contain  four  such  decifeet  components,  while  most  contain  three.  In 
figure  7 those  increments  with  four  components  have  been  modified  (with 
dashed  lines)  so  that  their  probabilities  conform  with  the  increments 
which  have  only  three  components.  It  can  be  seen  that  the  2. 4-2. 5 m 
peak  was  really  an  artefact,  while  the  2.0-2. 1 m,  2. 7-2. 8 m and 
5. 1-5.2  m peaks  remain  as  valid  peaks  of  probability. 

(’an  we  say,  then,  that  2.0-2. 1 m,  2. 7-2. 8 in  mid  5. 1-5.2  m 
represent,  say,  the  drafts  reached  by  first-,  second-  mid  multi-year 
undeformed  ice?  According  to  present  thermodynamic  theories  of  ice 
growth  in  the  Arctic  Ocean,  we  cannot.  The  results  of  Mavkut  and 
Intersteiner  (1071),  as  modified  and  tabulated  by  Thorndike  et  al 
(1075),  show  that  ice  growing  from  open  water  at  the  end  of  summer 
will  reach  a thickness  of  1.7(>  m by  April  10  of  the  following  year, 
and  2.04,  2.21  and  2.55  m by  April  10  of  its  second,  third  and  fourth 
years  of  growth.  All  these  values,  mid  the  yearly  depth  increments 
between  them,  are  less  than  the  level  ice  values  that  we  have  found. 

The  identification  of  level  ice  draft  with  ice  type  therefore  remains 
an  open  question. 
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FIGURE  b.  Probability 

i 

density  function  of  level 
ice  draft  for  '50-km  sections. 
Pep tb  increment  10  cm. 
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Mean  draft 

% 

level  ice 

Th  in 
ice  peak 

Level  ice  drafts  with 
probability  density  >0,5 

1 

3.81 

33.0 

0.0 

•) 

3.  17 

’.3.3 

2.0b 

32.3 

1 .0 

4 
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#5 
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X 
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0.8 
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/ 
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0.8 
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8 
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2.0 
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17 
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10 
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20 
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23 
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24 

2.07 
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3.61 

47.0 

1.0 

07 

— i 

. 

i % 

40.4 
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Mean 

3.07 

34.3 
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* Caribou  crossings 


The  values  of  mean  level  ice  draft  ill  table  111  are  not  very 
meaningful , because  some  ice  from  keel  bottoms  appears  in  the 
histograms  of  figs.  t>  and  7 and  biases  the  mean  draft  upwards. 

However,  the  figure  for  percentage  of  level  ice  in  table  111  is 
strongly  indicative  of  the  real  percentage  of  undeformed  ice  present, 
although  necessarily  an  overestimate  because  of  keel  bottoms.  Aguin 
we  see  a low  figure  for  sections  1 and  2 with  their  heavy  ridging;  a 
fairly  steady  equilibrium  value  for  sections  5 to  2*5  (which  passes  a 
run  test  for  homogeneity);  and  a final  decrease  of  level  ice  occurrence 
in  sections  2b  and  27.  The  thin  ice  peak,  as  shown  in  fig.  7 and 
table  111,  is  concentrated  either  around  0.2  to  0. 4 m or  around 
0.8  to  1.1  m.  Fig.  7 shows  0.8-0. 9 m as  the  most  probable  category, 
and  the  distribution  of  probabilities  in  the  depth  range  0-1.5  m in 
fig.  7 >s  a measure  of  the  relative  frequencies  of  polynyas  of 
varying  ages. 


5.  DISTRIBUTION  OF  KEEL  SPACINGS 


5.1  Independent  keels 

The  extent  of  an  independent  keel  is  defined  using  the  criterion 
that  the  troughs  on  either  side  of  the  keel  crest  (point  of  maximum 
draft)  must  descend  at  least  half  way  towards  the  local  level  ice 
surface,  in  this  case  defined  arbitrarily  as  a draft  of  2.5  m.  This 
is  analogous  to  the  Rayleigh  criterion  for  resolving  spectral  lines  in 
optics  and  is  identical  to  that  used  by  Williams  et  «1  (1975)*  Wadhams 
(1976,  1977b)  and  Weeks  et  «*1  (1977)  for  the  analysis  of  submarine  and 
aircraft  profiles  (for  airborne  laser  profiles  the  identification  of 
the  "local  level  ice  surface"  is  much  easier).  It  differs  from  the 
criterion  of  Hibler  et  al^  (197*0*  where  the  troughs  must  descend  a 
fixed  distance  (61  cm  for  surface  ridges)  from  the  peak;  Hibler  (1975) 
has  discussed  the  effect  of  this  difference  in  definition  on  the 
resulting  distribution. 


5.2  Theory  of  spacings 

Hibler  et  al^  (1972)  showed  that  if  ridges  occur  at  random  along 
a track  the  distribution  of  spacings  between  ridges  is  given  by 

Pr(x)  dx  = nexp  (-/u)  dx  ( 
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whore  /*  is  the  mean  number  of  ridges  per  unit  length  of  track  and 
1*  (x)  dx  is  the  probability  that  a given  spacing  lies  between  x and 
( \ + dx)  in  length.  Mock  et_  a_j^  (197-)  tested  this  relationship  for 
surface  ridges  using  aerial  photographs,  and  found  good  agreement  except 
for  an  excess  of  ridges  at  small  spaeings.  On  a purely  random  theory, 
however,  we  expect  a deficit  of  ridges  or  keels  at  small  spaeings,  on 
account  of  the  so-ca  lied  "ridge  shadowing"  effect  (Wudhums,  1077b). 

This  occurs  because  keels  have  a finite  slope  angle  so  that  their  crests 
cannot  lie  closer  than  a certain  minimum  distance  x(,rjt.  Within  this 

distance  the  shallower  ridge  is  not  detected  and  the  ridge-picking 
criterion  selects  only  the  deeper  ridge.  Figure  8 illustrates  this  effect 
for  two  keels  of  relief  h,  h'  (h'^h)  relative  to  the  level  ice  bottom, 
each  ridge  being  of  triangular  cross-section  w l th  slope  ot.  Under  these 
circumstances 

A theoretical  treatment  of  the  modification  of  (2)  by  (9)  is  complex, 
but  an  iipproximate  solution  is  given  in  a paper  awaiting  publication 
(Lowrx  and  Wadhams,  unpub  1.).  The  net  effect  is  that  close  spaeings 
and  shallow  ridges  tend  to  be  lost  preferentially  from  the  distributions 
of  spacing  and  draft. 


c o t o. 


(3) 


a. 9 Resu  1 t s 


I'lie  complete  numerical  results  for  keel  spaeings  are  given  in 
Appendix  l’.  The  results  are  presented  for  each  of  the  9U  km  sections, 
with  ;t  spacing  increment  of  20  ni.  Distributions  are  presented  for  all 
keels  deeper  thiui  9 m and  for  all  keels  deeper  than  9 m.  This  is  because 
Wadhams  (1977b)  found  that  a number  of  deep  floe  bottoms  appeared  in  the 
draft  range  9-9  m and  that  the  theoretical  keel  draft  distribution 
function  was  valid  only  beyond  9 m.  It  was  felt,  therefore,  that  by 
taking  9 m as  a cutoff  a more  valid  distribution  of  spaeings  of  "real" 
keels  could  be  obtained.  The  results  from  Appendix  C were  added  together 
to  yield  an  overall  keel  spacing  distribution  for  the  whole  profile;  this 
is  shown  in  figure  9. 

Both  distributions  (?  9 m and  >9  m)  show  general  agreement  with  (2), 
with  the  expected  deficit  tit  small  spaeings.  This  deficit  shows  itself 
only  in  the  spacing  range  0 - 4U  m,  as  opposed  to  the  results  of  Wadhams 
( 1977b,  fig.  4),  where  the  deficit  extends  its  influence  to  120  m.  This 
is  probably  because  of  the  transducer  boamwidth  in  the  SOVERFI GN  profile, 
which  makes  ridges  seem  broader  and  less  steep  than  they  really  are 
(Wadhams,  1077c).  At  large  spaeings  fig.  9 shows  a positive  deviation 
from  (2),  which  must  bo  due  to  an  additional  effect  upsetting  the  purely 
random  distribution.  1 suggest  that  this  effect  is  simply  the  presence 
of  leads  and  polvnyas,  which  interpose  occasional  smooth  stretches  of  ice 
into  the  otherwise  random  icefield  and  thus  generate  an  anomalous  number 
of  large  keel  spaeings. 

It  should  be  noted  that  the  lines  of  best  fit  to  the  rectilinear 
parts  of  fig.  9 do  not  have  a gradient  of  . For  the  9 m cutoff  the 
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NUMBER  OF  OCCURRENCES 


FIGURE  9.  Distribution  of  keel  spacings  over  whole  submarine  track.  Bin 
size  20  m.  Results  are  plotted  for  keels  deeper  than  5 m and  9 m,  and  a 
straight  line  is  fitted  to  the  central  portion  of  each  curve. 
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gradient  is  ^.2  ( = 7 • ’ ) while  for  9 m it  is  2.9 

( I.?).  One  expects  u gradient  of  magnitude  greater  than  m. 

because  keel  shadowing  implies  that  the  original  "population"  has 
been  reduced  before  entering  the  statistics.  This  is  so  for  9 m 
but  not  for  *>  in,  again  implying  that  9 in  is  a better  cutoff  to  use 
so  a>  to  obtain  a population  of  "pure"  keels. 

Appendix  (’  also  presents  a tabulation  of  spacings  classified 
according  to  the  depth  of  the  deeper  keel  of  the  pair  which  defines 
the  spacing,  i.e.  a tabulation  of  x against  (h1  + 2.‘i)  in  figure  8. 

This  is  an  attempt  to  find  the  best  cx  for  use  in  (A).  In  Appendix  C 
the  tabulation  is  given  for  all  27  sections  and  for  A in  and  9 m 
cutoff.  We  have  taken  the  results  for  9 m cutoff  and  plotted  them 
in  bins  of  2 m depth  increment  and  20  m spacing  increment.  The  curves 
(figure  111)  show  a peak  at  a sparing  which  progressively  increases  with 
depth.  In  fig.  10  the  position  of  the  highest  peak  in  each  curve  has 
been  plotted  against  the  relevant  depth  f the  heavy  black  dots),  and  it 
can  be  seen  that  the  increase  with  depth  is  roughly  linear.  A line  of 
best  I’ll  has  been  drawn  through  these  points  which,  when  applied  to  (A), 
gives  a value  of  11. V for <x . Of  course,  these  curves  show  that  there 
is  no  one  value  of  cx  , otherwise  there  would  be  a sharp  spacing  cutoff 
within  which  no  keel  pairs  are  to  be  found.  Instead,  there  is  a range 
of  <x,  a range  which  is  spread  out  still  further  by  the  fact  that  the 
keels  are  not  being  profiled  orthogonally  but  at  various  angles  of 
encounter.  Wadliams  (1977c)  dealt  with  this  statistical  averaging 
problem.  Our  value  of  11.1°  is,  in  any  case,  an  underestimate 
because  it  refers  to  the  peak  of  each  spacing  distribution  rather  than 
to  the  spacing  at  which  keel  pairs  begin  to  be  found.  However  it  is 
indicative  of  the  validity  of  the  "keel  shadowing"  concept.  Again, 
Wadhams  (1977c)  found  by  actually  measuring  the  slopes  of  keels  on  sonar 
profiles  that  the  slope  angle  distribution  had  a peak  in  the  range 
10-20°  but,  when  adjusted  to  take  account  of  imgle  of  encounter,  the 
mean  value  of  <x  came  to  A2°.  Thus  our  1A.A°  figure  does  not  appear 
unreasonably  low. 


0.  D1STH1B11  l lON  OF  KEK1.  DRAFTS 

0. 1 Theory 

The  theory  of  keel  drafts  which  has  been  most  extensively  tested 
against  observation  is  that  of  Hibler  e_t_  iU  (1972).  They  used  a 
variational  calculation  which  gives  the  most  likely  distribution  of 
geometrically  congruent  ridges  that  will  yield  a given  volume  of 
deformed  ice.  The  result  is 

I’r(h)  dh  = 2Ah  exp  (Mijj)  exp(->h“)  dh  (4) 

where  I’r(h)  dh  is  the  probability  that  the  draft  lies  between 
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SPACING  m 


FIGURE  10.  Spacing  distributions  for  keels  of  draft  greater  than  9 m, 
in  2 m draft  increments.  The  spacing  corresponding  to  the  maximum  of 
each  distribution  has  been  plotted  as  a heavy  black  dot.  This  spacing 
increases  linearly  with  keel  relief. 
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h and  (h  + iili ) ; 
h is  the  mean  draft; 

h()  is  a low  value  cutoff  below  which  keels  are  not 
included  in  the  statistics; 

\ is  a parameter  which  must  be  derived  by  iteration  from 

exp(->h“)  = h ( XT!) - erfc(>-  h0)  (5) 

This  has  been  shown  to  give  a better  fit  to  submarine  sonar 
observations  than  an  alternative  distribution  proposed  on  empirical 
grounds  by  Diachok  (1975): 

Pr(h)  dh  = ~tt  exp  (-h“/a^)  dli  (b) 

tl" 

with  a = 2h/n- . 

Again,  (’i)  is  modified  by  the  ridge  shadowing  effect  (Lowry  and 
Wadlioms,  unpub  1.),  but  the  modification  has  only  a small  effect  at  the 
low-draft  end  of  the  distribution  where  it  causes  a slight  deficit  of 
keels;  this  is  a much  less  drastic  modification  than  that  applied  to 
the  spacing  distribution. 

Recently  it  has  been  found  that  surface  ridge  sails,  to  which  this 
theory  was  also  thought  to  apply,  actually  obey  a simpler  negative 
exponential  distribution  of  form 

l’r(h)  dh  = B exp  (-bh)  dh  (?) 

with  B,b  as  parameters,  provided  the  sails  are  identified  using  the 
Rayleigh  criterion  (Wadhams,  197b;  Weeks  et  al , 1977)*  Hi  bier  (1975) 
showed  that  the  same  data  can  be  made  to  fit  or  (7)  depending  on 

whether  the  Hibler  (constant  trough  depth)  or  Rayleigh  ridge-picking 
criterion  is  used. 


b.2  Results 


Figure  11  shows  the  distribution  of  keel  drafts  for  all  27 
sections  at  50  km  gauge,  expressed  ns  keels  per  100  km  track  and  using 
a 1 m depth  increment.  The  first  bin  (2-3  m)  actually  contains  only 
keels  from  2.5  m to  3 m in  draft,  since  2.5  m is  the  zero  datum  for  the 
Rayleigh  criterion.  "Keels"  of  less  than  5 m draft  can  be  assumed  to 
consist  mainly  of  the  bottoms  of  undulating  floes. 

The  overall  data  from  all  27  sections  were  collected  and  are 
plotted  in  figure  12  on  a semi-log  scale  for  all  keels  deeper  than 
5 m.  The  result  is  an  exceptionally  good  fit  to  a straight  line,  not 
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FIGURE  it.  Distribution  of 
keel  drafts  for  50-km  sections. 
Depth  increment  1 m. 
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FIGURE  12.  Distribution  of  keel  drafts  plotted  on  a semi-log 
scale  for  data  at  1*400  km,  400  km  and  *50  km  length  gauges. 
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only  for  the  overall  data  hut  also  for  data  at  1200  km  gauge  (e.g. 
section  F,  also  shown  in  fig.  Ill)  tuul  even  'll)  km  gauge  (section  1, 
fig.  12).  This  shows  that  the  keel  draft  distribution  obeys  the 
simpler  relationship  (7)  rather  than  the  Ilibler  relationship  (4). 

This  is  a most  unexpected  result,  because  (?)  was  hitherto  thought 
to  be  valid  only  for  ridge  sails;  the  SOVEREIGN  keel  data,  analysed 
using  the  Rayleigh  criterion,  follow  (4)  with  u high  degree  of 
exactness  as  do  all  other  published  sonar  profiles.  Wudhams  ( 1077b) 
suggested  that  sails  may  not  follow  (4)  because  they  contain  only  a 
small  proportion  of  the  mass  of  a ridge  and  their  shape  is  determined 
largely  by  accident.  This  cannot  explain  the  present  result.  We 
must  conclude  either  that  keels  in  the  Beaufort  Sea  have  a different 
nature  from  those  in  the  Eurasian  llasin,  which  is  unlikely,  or  that 
the  apparent  distribution  of  independent  keel  drafts  is  dependent  in 
some  way  on  the  type  of  sensor  employed. 

Profiles  of  ridge  sails  are  always  obtained  using  a laser 
prof i lome ter , which  has  a pencil  beam  capable  of  recording  much  of  the 
fine  structure  of  the  sail,  including  crevices  and  troughs  between  the 
blocks  (the  limitation  being  the  integration  time  of  the  laser  electronics). 
The  sonar  employed  by  GURNARD  also  had  a narrow  heamwidth  (its  precise 
characteristics  being  classified)  and,  as  shown  in  fig.  T,  it  also  appears 
capable  of  recording  the  fine  structure  of  a keel,  probing  into  clefts 
and  hollows  between  the  submerged  blocks.  DREADNOUGHT  (Williams  e_t_  al , 

1 1)7‘5 ) used  a sounder  with  a very  wide  bo  tun , and  SOVEREIGN  (Wailhams, 

1977b)  had  a sounder  with  a wide  beam  in  the  fore-and-aft  plane  (17°) 
and  a narrow  beam  in  the  athwartships  plane  (*>  ).  Wiile-beam  sounders 
smooth  out  the  structure  of  a keel  so  that  it  is  always  perceived  as  a 
single  wedge  (see,  for  example,  the  profiles  in  Wadhtuns,  1077c),  and 
even  the  application  of  reconstruction  equations  (Williams  ej_  nl_,  197^) 
cannot  regenerate  this  fine  structure.  Now  Ilibler' s theory  depends  on 
the  concept  of  geometrically  congruent  ridges,  each  iui  entity  of  the 
same  shape  possessing  mass  and  potential  energy  which  depend  only  on  its 
depth.  A wide-beam  sounder  forces  keels  to  approximate  to  this 
concept  by  smoothing  out  any  incidental  structure  that  they  may  possess 
and  leaving  them  as  discrete  entities.  Thus  narrow-beam  echo  sounders 
and  laser  prof i lometers  produce  one  type  of  ice  profile  with  ridge 
height  characteristics  obeying  (7)*  while  wide-beam  sounders  produce 
another  type,  obeying  (4).  A narrow-beam  sounder,  by  splitting  many 
ridges  into  multiple  "ridges",  sees  a greater  ridge  frequency  than  a 
wide-beam  sounder  (e.g.  in  Wailhams,  1977b,  the  sail  frequency  is  a 
multiple  of  the  keel  frequency  for  the  same  ice  cover). 

This  hypothesis  now  covers  all  results  except  laser  prof i lometer 
data  reported  by  Ilibler  e_t_  a_l  (1974),  which  still  obeyed  (4)  but  which 
were  not  analysed  on  the  Rayleign  criterion.  A crucial  test  of  the 
hypothesis  would  be  to  smooth  the  GURNARD  profile  artificially  by 
convolving  it  with  the  beam  pattern  of  a wide-beam  echo  sounder,  anil 
to  observe  the  effect  on  the  resulting  statistics.  We  hope  to  report 
on  this  computer  simulation  in  a later  publication. 

It  was  shown  in  Wadhams  (1977b)  that  the  parameters  B,b  in  (7) 
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cun  be  expressed  as  simple  functions  of  h and  /a,.  If  (7)  is  rewritten 
in  the  form 


n(h)  dh  = B exp  (-b  h)  dh 

where  n(h)  is  the  number  of  keels  per  km  track  per  metre  draft 
increment,  then 


and 


M.  = J n(h)  dh 
ho 


H = Jin 


(h)  dh 


(8) 


(9) 


(10) 


so  that 


and 


b = (h  - h J 


-1 


B = ycb  exp(b  h0) 


(11) 

(12) 


Thus h and  h are  the  two  parameters  of  the  keel  draft  distribution 
from  which  the  whole  shape  of  the  distribution  can  be  deduced  using  (8) 
to  (12).  Table  IV  shows  these  parameters  tabulated  for  all  the  50  km 
sections,  using  h = 5 m and  9 m so  as  to  be  consistent  with  the 
statistics  of  Wadhams  ( 1977b).  The  results  are  also  plotted  in 
figure  15.  The  following  conclusions  can  be  drawn:- 

(a)  The  major  part  of  the  profile  (sections  5-25)  has  a ridging 

distribution  which  is  very  homogeneous  and  which  falls  within 
narrow  limits  of  variation.  These  limits  are  extremely 
narrow  for  h0  = 5 m (5. 4-8. 2 for/*.;  7*2-7«8  m for  h)  and 

somewhat  wider  for  h0  = 9 m,  probably  because  of  the  smaller 
number  of  keels  involved.  Mean  and  standard  deviation  for 
these  four  parameters  are  given;  run  tests  show  that  we  can 
accept  the  hypothesis  of  a homogeneous  ice  cover  with  respect 
to  ridging  intensity. 

(b)  Sections  1-2  and  26-27  fall  clearly  outside  the  range  of 
variation  of  the  other  sections  (fig.  15),  indicating  much 
heavier  ridging  (greater/*. ) in  these  parts  of  the  track. 

This  result  agrees  with  what  we  have  found  from  the  probability 
density  functions;  at  these  two  extremities  of  the  track  there 
are  more  pressure  ridges  per  unit  length,  a greater  mean  keel 
draft,  a greater  mean  ice  draft  and  a greater  proportion  of 
deformed  ice. 

(c)  There  is  no  clear  positive  correlation  between  /*.  and  h, 
although  the  four  heavily  ridged  sections  have  both  a high  /*. 
and  a high  K. 
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Table  IV. 

Pressure  ridge 

frequent’  l es 

and  mean 

draf  t s 

for  50  km 

sect  ions 

Sect i on 

Draft  > 1 m 

No.  per  km  Mean  draft 

Draft  >9 
No.  per  km  Mt 

m 

•an  draft 

Maximum 
draft  m 

1 

i i . < < 

8.57 

5.  1 1 

12.15 

21.11 

•» 

10.28 

8.10 

1.07 

1 2 . 20 

28.81 

8.  1 1 

7.t>o 

1.81 

1 1 . 7‘> 

21.04 

4 

7.41 

7.02 

1.55 

l 1 .75 

19.00 

5 

7.02 

7.78 

1 .00 

1 1 . 0 1 

22.05 

0 

1 . t>8 

i • 5 / 

1 . 10 

1 •>  77 

22.05 

t 

1 . 49 

7 .41 

0 . 00 

12.00 

20.54 

8 

1 .'*  7 

7.20 

0.01 

1 1 .84 

20.73 

0 

8.  17 

7.10 

1 .41 

11.51 

20.00 

11) 

0.40 

7.51 

1 . 40 

1 1 .14 

10.80 

1 1 

0.50 

7.2t) 

1 . 24 

1 l . 25 

22.59 

12 

7-08 

7-b5 

1.58 

11.75 

23.05 

(Short 

t>.  17 

7.78 

1 .00 

11.  10 

11.12) 

n 

ti . bO 

7.b5 

1.40 

11.80 

20.07 

14 

b.bO 

7.08 

1.18 

12.  11 

24 . 91 

15 

5.4b 

7.12 

0.80 

12. 1 1 

21.18 

lb 

0.24 

7.19 

1.11 

11.23 

18.90 

17 

7.15 

7.7‘> 

1.05 

12.11 

24.84 

18 

7.4  b 

7 . 04 

1 . 04 

1 1 . 40 

18.15 

10 

1 . 10 

7.51 

1.11 

1 1 . 91 

20.71 

20 

7 or 

7.50 

1.18 

11.71 

20.10 

21 

b.98 

7.52 

l . 51 

11.18 

19.99 

')•> 

5.41 

/ * i 1 

1.17 

11.71 

19.00 

21 

7.01 

7.05 

1.50 

12.04 

24 . 48 

2’* 

b.lb 

r . 40 

1.22 

12.01 

22.01 

25 

7.62 

7 . 00 

1.78 

11.72 

22.07 

2b 

8.82 

8.25 

2.05 

12.12 

20.21 

07 
— 1 

8.72 

8.14 

2.00 

12.09 

29. 14 

Mean  of  main 

6.71 

1 .10 

11.79 

liroup  (1-25) 

+0.20 

+0.03 

+0.00 

+0.07 

Mean  of  Caribou 

crossings 

b.bl 

7.00 

1.42 

11.88 

( *»— 7 * 20-21) 

SOVEREIGN 

5.08 

11.74 

1.74 

14.10 

41+1 

DREADNOUGHT 

4 . 20 

9.57 

2.00 

12.57 

10+ 1 
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FIGURE  13.  Mean  keel 
draft  plotted  against 
mean  number  of  keels 
per  km  for  keels  deeper 
than  3 m and  9 m.  Each 
point  represents  a 
30  km  section. 
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(d)  For  h„  = l)  m both  the  frequency  the  mean  draft  are  much  lower 

than  those  found  b\  SOVEREIGN  in  the  very  heavily  ridged  zone  off 
north  Greenland.  Only  section  1 exceeds  the  SOVEREIGN  data  in 
keel  frequency,  though  not  in  mean  draft.  The  data  compare  very 
well  with  DREADNOlHiHT  data  from  the  central  Eurasian  Basin: 
although  it  may  appear  from  table  IV  that  DRKADNOUtiHT  data  have  a 
lower /<  and  higher  h,  the  difference  can  be  ascribed  to  the  progrum 
of  Will  nuns  ej^  a_l  (IO73)  which,  by  applying  a harsh  version  of  the 
Rayleigh  criterion  with  sea  level  as  the  zero  datum,  lost  many 
shallow  keels  from  the  statistics. 


1 
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(e)  The  maximum  drafts  are  surprisingly  low.  On  the  whole  the  deepest 
keel  drafts  are  found  in  the  four  anomalous  sections  (1-2,  26-27), 
but  the  deepest  draft  of  all,  31.12  m,  in  fact  occurred  in  the  short 
18  km  portion  of  track  that  was  omitted  from  the  30  km  statistics. 

The  general  ridging  properties  of  this  portion  (tablelV)  are  quite 
typical  of  the  sections  surrounding  it,  so  that  the  keel  can  be  seen 
as  an  isolated  event.  This  is  the  only  keel  deeper  than  30  m in  the 
entire  1400  km  of  profile,  whereas  in  the  SOVEREIGN  profile  there  were 
43  keels  deeper  tlwui  30  m in  3000  km  of  track,  30  of  them  occurring  in 
the  1030  km  of  "offshore  zone"  north  of  Greenland  (Wadhams,  1077c). 

By  a coincidence  the  deepest  keel  in  the  SOVEREIGN  profile,  43  m,  also 
occurred  as  an  isolated  event  in  an  otherwise  lightly  ridged  section 
of  ice  cover. 
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LEADS  AND  P0LYNYAS 


The  probability  density  function  of  ice  draft  gives  the  best  measure 
of  the  occurrence  find  thickness  distribution  of  the  thin  ice  in  leads  and 
polynyas,  and  is  especially  useful  for  application  to  heat  budget  calcula- 
tions. However  it  is  important  for  a variety  of  applications  in  ice 
mechanics,  traff icabi 1 ity  etc.  to  know  the  frequency  and  width  distribution 
of  leads  encountered  by  the  submarine.  Perhaps  the  most  important 
application  is  to  submarine  operations  themselves  - it  is  desirable  to  know 
the  mean  spacing  of  leads  that  are  large  enough  to  permit  a submarine  to 
surface. 


A lend  was  defined  as  a continuous  sequence  of  depth  points  in  which 
no  point  exceeds  1 m in  draft  - thus  a polynya  broken  up  by  a small  floe 
of  broken  ice  counts  as  two  leads.  Lead  widths  were  classified  in  30  m 
increments  and  the  results  calculated  for  each  of  the  three  files  making 
up  the  overall  track  (File  l = sections  1-12;  File  2 = 13-21;  File  3 = 
22-27).  The  results  are  shown  in  table  V. 

On  average  the  aquatic  crow  has  to  swim  only  about  200  m between  leads  - 
although  this  figure  varies  by  a factor  of  nearly  4 between  File  2 find 
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Table  V (a) 

• 

i 

> 

Distribution  oj 
encountered 

lead  widths 
per  100  km  of 

: number  of 

track. 

1 eads 

Lead  width  in 

ri 

File  1 

File  2 

File  3 

Overa 

0-30 

389 

257 

915 

453 

30-100 

5.6 

5.2 

4.4 

5.6 

100-150 

0.9 

1.5 

2.8 

1.5 

150-200 

0.5 

2.1 

0.9 

1.1 

; 200-250 

0.3 

0.2 

0.3 

0.3 

250-300 

0.6 

0.5 

0.5 

0.6 

300-350 

0 

0 

0 

0 

350-400 

0.2 

0 

0 

0.1 

400-450 

0 

0 

0 

0 

450-500 

0 

0.2 

0 

0.1 

500 

0.5 

0.4 

0 

0.4 

I 
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(b)  Mean  distance  travelled  between  leads 


UidtL  m Distance  between  encounters 

0 212  m 

‘50  10.3  km 

100  24.1  km 

130  38 • 4 km 

200  b?.6  km 

300  237  km 
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File  5.  However,  the  lead  which  he  reaches  is  likely  to  be  very 
narrow,  and  few  leads  exceeded  50  m in  width.  The  "exceedence 
table"  in  V(b)  shows  that  a submarine  which  requires  a 200  m lead 
for  a safe  surfacing  will  have  to  travel  b8  km  to  find  one.  In 
fact  a submarine  trying  to  surface  usually  investigates  every  lead 
wider  than  about  50  m in  case  the  cross-track  dimension  is 
sufficient  to  permit  surfacing;  50  m leads  occurred  every  10  km 
in  the  southern  Beaufort  Sea  during  the  period  of  this  experiment. 


APPENDIX  A.  ARCTIC  SUBMARINE  LABORATORY  PROCESSING  TECHNIQUE 

The  digitized  magnetic  tapes  were  unpacked  and  processed  by  the 
Arctic  Submarine  Laboratory  on  a ITNIVAC  1108  - 1110  computer  system 
at  the  Naval  Undersea  Center.  Time  marks  were  removed,  and  the 
corrected  data  was  then  merged  with  position  information  and  converted 
to  BCD  format  for  transmittal  to  Project  AIDJEX. 

After  the  initial  unpacking,  spurious  profile  points  resulting  from 
multiple  echoes,  fish,  air  bubbles,  etc.,  were  eliminated  by  two  correction 
programs.  Any  points  that  were  deleted  were  set  to  -10.0.  Multiple 
adjacent  deleted  points  were  left  at  this  value,  but  single  points  were 
filled  by  linear  interpolation.  These  were  flagged  by  adding  200  to  the 
interpolated  value. 

The  first  correction  program  deleted  any  points  that  were  greater 
than  or  equal  to  150  feet.  At  the  beginning  of  a tape,  or  following 
one  or  more  zeros  recorded  by  the  system,  the  program  searched  for  the 
first  valid  data  point  by  computing  the  absolute  difference  between 
adjacent  points  (|X.  - X^+^|),  and  requiring  that  the  difference  be  less 
than  5.0  feet.  For  example,  if  the  difference  between  the  first  two 
points  satisfied  this  criterion,  the  program  assumed  that  X^  was  correct 
and  used  it  as  its  current  X.  value.  If,  however,  the  difference  was 

l 

greater  than  5.0,  it  deleted  X^,  made  the  current  X^  value,  and 
computed  a difference  (|  Xg  - X^l).  It  continued  in  this  manner  until  it 
found  a difference  that  was  less  than  5.0,  at  which  time  it  accepted  the 
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current  X.  value  as  the  first  valid  data  point. 

After  the  program  chose  a valid  X.  value,  it  conqmted  the  absolute 

difference  between  adjacent  points  luid  compared  the  difference  with  a 

threshold  number  that  was  set  to  20.0  feet.  If  the  difference 

(JX^  - X.+j|)  was  less  than  the  threshold,  the  program  accepted  X^  + j 

as  a valid  point  and  used  it  as  its  next  X.  value.  However,  if  the 

difference  was  greater  than  20,  it  computed  a I),,  equal  to  |X^  - X^Q|, 

and,  if  necessary,  a 1)^  equal  to  |X^  - X.  _|  . The  program  uccepted 

X.  + | as  a valid  point  only  if  both  D,,  and  1)^  were  also  greater  than  the 

threshold  value.  If  I),,  was  less  than  20,  X.  . was  deleted,  and  X.  „ 

was  selected  as  the  next  X.  value.  If  1)1}  was  greater  than  20  but 

1)_  was  less,  both  X.  . mid  X.  ,,  were  deleted,  mid  X.  _ was  chosen  as 
1 1+1  i+2  i+i 

the  next  X.  value. 

i 

Since  the  system  recorded  zeros  during  periods  when  it  didn't 

receive  any  data,  zero  words  were  not  used  to  compute  the  differences. 

The  program  required  that  if  I) ^ ( | X ^ was  greater  than  the 

threshold,  both  D1}  mid  must,  be  greater  tliun  20  or  the  data  were 

rejected.  In  the  case  that  D ^ was  greater  than  20  and  X.+1)  was  zero, 

the  program  deleted  X.  + j If  Dj  mid  0,,  were  greater  than  20,  but. 

X.  „ was  zero,  both  X.  . and  X.  were  deleted.  In  either  case,  or  if 
i+>  i+l  i+2 

X . + 1 was  zero,  the  program  searched  for  the  next  valid  X.  value  by  using 
the  technique  described  ubove  (third  paragraph). 

In  order  to  maintain  continuity  between  duta  records,  the  program 
included  the  first  three  data  points  from  the  succeeding  record  with 
the  current  record  when  it.  computed  the  differences. 

The  second  program  subtracted  a surface  offset  correction  from 
each  of  the  profile  points.  This  correction  was  determined  from  mi 
analysis  of  both  analog  mid  digital  data,  mid  is  used  to  shift  the  mean 
data  points  in  areas  of  open  water  to  read  zero.  A few  points  that,  were 
already  very  near  zero  became  negative.  The  offset  was  usually  one  to 
two  feet. 
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This  program  also  eliminated  .spurious  points  that  might  have 
been  missed  by  the  first  program  by  deleting  any  points  during  a 
variable  time  period  that  exceeded  the  maximum  possible  ice  depth 
which  wus  determined  by  examining  the  unulog  records.  Points 
deleted  bv  this  program  were  not  interpolated  and  appear  us  -10.0. 


APPENDIX  II.  Si'll  1 COMPUTER  PROGRAM  FOR  DATA  PROCESSING 


The  program  is  written  in  Algol  68C.  It  consists  of  three 
parts:- 


(a)  the  main  program,  which  controls  the  whole  computation, 
reading  in  data,  accumulating  statistics  and  outputting 
results; 

(b)  a series  of  statistics  procedures  for  initialising,  updating 
and  outputting  tables,  identifying  level  ice  and  ridges  etc. 
Note  that  two  duta-structures  are  defined,  MODE  TABLE  and 
MODE  TWOWAYTABLK.  Each  of  these  structures  contains  details 
about  the  classification  defining  the  table,  as  well  as  space 
for  storing  the  accumulated  frequencies  etc.; 

(c)  various  miscellaneous  "utility"  procedures,  which  may  be 
considered  as  just  an  extension  to  the  Algol  68C  stnndurd 
prelude. 
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APPENDIX  C.  HILI,  NUMERICAL  RESl’LTS 


The  numerical  results  from  the  foregoing  program  for 
all  27  sections  are  lengthy  and  have  therefore  been  published 
separately  as  a Supplement  to  this  Report.  Readers  who  have 
not  received  a copy  of  this  Supplement  with  their  Report  may 
obtain  one  by  writing  to  the  authors  at  Scott  Polar  Research 
Institute,  Cambridge  CB2  1ER,  Great  Britain. 
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